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While large scale cosmic microwave background (CMB) anisotropics involve a combination of the 
scalar and tensor fluctuations, the scalar amplitude can be independently determined through the 
CMB-galaxy cross-correlation. Using recently measured cross-correlation amplitudes, arising from 
the cross-correlation between galaxies and the Integrated Sachs Wolfe effect in CMB anisotropics, 
we obtain a constraint r < 0.5 at 68 % confidence level on the tensor-to-scalar fluctuation amplitude 
ratio. The data also allow us to exclude gravity waves at a level of a few percent, relative to the den- 
sity field, in a low - Lambda dominated universe (Oa ~ 0.5). In future, joining cross-correlation ISW 
measurements, which captures cosmological parameter information, with independent determina- 
tions of the matter density and CMB anisotropy power spectrum, may constrain the tensor-to-scalar 
ratio to a level above 0.05. This value is the ultimate limit on tensor-to-scalar ratio from temper- 
ature anisotropy maps when all other cosmological parameters except for the tensor amplitude are 
known and the combination with CMB-galaxy correlation allows this limit to be reached easily by 
accounting for degeneracies in certain cosmological parameters. 

PACS numbers: 



I. INTRODUCTION 

One of the major predictions of inflation is the exis- 
tence of a stochastic background of gravitational waves 
(GW) (see e.g. pj, Q)- The amplitude of these tensor 
fluctuations is proportional to the square of the energy 
scale of inflation (see e.g. 3]). Furthermore, the 'tilt' 
of the GW spectrum (and of the scalar as well) can give 
us direct information on the slope of the inflaton poten- 
tial up to second derivatives. Therefore the detection of 
this stochastic background would provide valuable infor- 
mation on the inflationary scenario shcding light on the 
physics at ~ I0 16 GeV (see e.g. 0). 

The stochastic GW background leaves a signature on 
the CMB through different mechanisms. In principle a 
clear detection can be obtained from measurements of the 
B-polarization. Tensor perturbations also induce tem- 
perature fluctuations through a Sachs- Wolfe effect and 
an integrated one (for a review see The former is 

localized by the visibility function at the last scattering 
surface and consequently depends on the recombination 
history. Since it occurs in a relatively short period of 
time the amplitude of this effect is negligible at all scales. 
On the contrary GW produce most of the signal through 
the integrated effect which is sensitive to the evolution 
of gravitational waves from the time of last scattering to 
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today. Since the observed temperature anisotropy power 
spectrum is the sum of both the scalar and tensor compo- 
nents, for relatively low amplitudes of the tensor modes 
the scalar fluctuations dominate the spectrum and the 
tensor contribution is dwarfed. Similarly, if there is a 
non-negligible contribution from GWs, the sum of scalar 
and tensor spectra would lower the predicted amplitude 
of the acoustic peaks relative to the angular scales of few 
or more degrees. Thus, the ratio of tens of degree to sub- 
degree fluctuations can be used to constrain the presence 
of GW. 

With the advent of WMAP, several studies have ad- 
dressed whether there is a significant contribution of 
GWs in the CMB spectra (see e.g. 0,0,1113). Despite 
the high quality of these data, it still remains difficult to 
obtain a robust constrain on the amplitude of the tensor 
perturbations due to the degeneracies with other cosmo- 
logical parameters. For instance, both the scalar spectral 
index ns and its running are strongly degenerate with the 
amplitude of the GW background. In [lfj , combining the 
CMB data with the 2dF matter power spectrum yield a 
2a constraint on the tensor — to — scalar ratio r < 0.9 
[39|. The combination of WMAP data with Sloan Ly-cv 
forest, lensing bias, and galaxy clustering data constrain 
r < 0.36 at the 95% confidence level in a flat cosmological 
model • Measurements of B-mode CMB polarization 
from future experiments such as the Planck satellite, or 
in the long-term the CMBpol mission as part of NASA's 
Beyond Einstein Program, will provide either a detec- 
tion of the tensor modes or the best constrain on their 
amplitude 12]. Concept studies are already under way 
to consider the possibility for a direct detection of the 
relic gravitational wave background present today with 
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space-based Laser interferometers as part of a second- 
generation mission to followup NASA's Laser Interfer- 
ometer Space Antenna (LISA; [13 )• In the meantime, it 
is worthwile to investigate methods that could indepen- 
dently bound the GW contribution and/or break some of 
the degeneracies which affect present constraints. This 
would also be useful in planning observational strategies 
of current and upcoming CMB experiments. 

The scalar and tensor temperature anisotropy power 
spectra are characterized by a nearly constant plateau 
on the large angular scales. The temperature anisotropy 
power specrum are sums of tensor and scale contribu- 
tions, Cj T = cf T ' scalar + cf T ' tensor . To establish the 
amplitude of tensor fluctuations from total temperature 
anisotropies, which is what we observe, one must need 
to properly account for the contribution from scalar fluc- 
tuations. Thus, an independent estimate of the scalar 
amplitude, when combined with temperature anisotropy 
spectrum, can be used to constrain the amplitude of ten- 
sor fluctuations, or, the ratio of tensor-to-scalar fluctu- 
ation amplitude. We should always keep in mind that 
due to the integrated Sachs- Wolfe effect produced at late 
time by the scalar perturbations, (jJ" T ' scalar has a non- 
trivial normalization that depends both on the potential 
at last scattering (the Sachs- Wolfe effect), as well as the 
time evolution of the potential; The normalization of the 
CMB temperature power spectrum also depends on the 
amount as well as the nature of the dark enegy given 
that potential fluctuations evolve differently depending 
on late-time cosmology 0, 0, . 

The scalar amplitude can be inferred from the matter 
power spectrum as measured by galaxy surveys [TtI Il8| . 
Alternatively we can make use of the positive cross- 
correlation between WMAP anisotropy maps and LSS 
surveys H3M 110. Here we consider such a pos- 
sibility and determine an independent constrain on the 
tcnsor-to-scalar ratio. Our limit improves those inferred 
from the combination of the WMAP data with either 
2dF or SDSS surveys 0,0, by nearly a factor of two. 
In particular we find that in a flat cosmological constant 
dominated universe r < 0.59 at 2a. We also discuss how 
the constraints on r change by including more accurate 
information on the ISW effect. 



II. LIMIT ON TENSOR-TO-SCALAR RATIO 

The correlation between the large angular scale CMB 
temperature anisotropy maps and the distribution of the 
large scale structure (LSS) is consequence of the inte- 
grated Sachs-Wolfe. As the universe starts deviating 
from the matter dominated expansion, the gravitational 
potentials associated with clumps of matter change in 
time. CMB photons crossing these regions undergo a 
shift which causes a temperature anisotropy. Since the 
clustered structures arises from the scalar fluctuations 
the CMB-LSS correlation depends only on the scalar 
modes and does not get any contribution from the tem- 



perature anisotropies induced by the gravitational waves 
background. Indeed the latter is uncorrelated with the 
distribution of galaxies, cluster and super-cluster. Thus, 
the cross-correlation between CMB maps and surveys can 
be used to infer A s , under the assumption of a given dark 
energy model. 

The angular cross-correlation in Legendre series is 
given by 

00 O/ll 

C x (9)=Y.^C?P l (co S (e), (1) 

1=2 

where Pi(cos9) are the Legendre polynomials and Cf- is 
the cross-correlation power spectrum 0] : 

C? = IJ k 2 dk P s& {k)ll sw {k)lt ss {k) . (2) 

Here, P$s(k) is the power spectrum of density fluctua- 
tions (which captures the amplitude of scalar fluctua- 
tions) and integrand functions lf sw (fc) and I^ ss (k) are 
defined as 

l{ SW (k) = J W lsw {k,r)ji[kr]dr (3) 
lt SS (k) = J W^ ss {k,r)ji[kr]dr (4) 
when the two window functions are 
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W LSS (k,r) = b(k,r)n{r)G{r). (5) 

In above, ji[kr] is the spherical Bcsscl function when r(z) 
is the comoving distance, 6(fc, r) is the scale- and redshift, 
or distance-dependent bias factor of galaxies with a nor- 
malized distribution, in distance, of n(r) (f drn{r) = 1), 
and G(r) is the growth function of density fluctuations. 
Other parameters have standard definitions. In convert- 
ing time evolving potential fluctuations, captured by the 
ISW effect, to density Auctions, we have assume a flat- 
comsological model. If not, the window function re- 
lated to the ISW term, W lsw (k,r), will contain cor- 
rections related to curvature. We refer the reader to 
Ref. [H 0, for a more detailed discussion 
on the cross-correlation power spectrum. 

From Eq. J3J, it is evident that the determination of 
A s from cross-correlation measurements depends on the 
galaxy bias b(k, z) of the LSS catalog used. Nevertheless, 
we can state in a conservative way that the r.m.s. value 
of the matter fluctuation at R = 8Mpc//i lies in the range 
<Tg = 0.7 — 1.1 15]. Normalizing the predicted matter- 
matter power spectra to these values and comparing to 
the observed galaxy-galaxy correlation function one can 
determine the corresponding bias. Marginalizing over as 
introduces, for a single catalog, an extra ~ 20% error in 
the measurement of the ISW-correlation signal. 
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FIG. 1: Likelihood contours in the Qa — r plane from current 
ISW detections (sec 33]). Flatness is assumed. The ISW sig- 
nal is smaller for small values of the Cosmological Constant. 



This error can be further reduced if one considers sev- 
eral ISW detections with catalogs affected by different 
biases, or, in the case of a single survey, if biasing can be 
determined as a function of the galaxy luminosity |29j . 
Here is worth remarking that the dependence of the cross- 
correlation on the bias is linear. Therefore, the error on 
the inference of A s induced by the uncertainty on the 
bias is fractionally less here than in the standard method 
involving the combination of CMB anisotropy data and 
LSS clustering spectrum (e.g. [13 )> as the full galaxy- 
galaxy correlation function quadratically depends on the 
galaxy bias factor. 

Another caveat concerns the dark energy. Current 
measurements of the cross-correlation lack the accu racy 
to strongly constrain any dark energy parameters |27| . 
The next generation of large scale structure surveys will 
provide better constraints by measuring the correlation 
between the CMB and different redshift sample of galax- 
ies [s^. As shown in [2^, [2^ the redshift behavior of 
the cross-correlation signal depends on the specific of the 
dark energy, while A s is an overall normalization factor 
independent of the redshift. 

The cosmological constant is the simplest dark energy 
model to best fit current data [Tol Il7l 111 l3ll E2} . 

We thefore limits our analysis to flat LCDM cosmolo- 
gies and use the current ISW detections as collected 
by [3^ . These data consist of average angular cross- 
correlation measurements between WMAP temperature 
maps and 5 different LSS surveys. For a more detailed 
discussion on the evaluation of the likelihood statistics, 
the effect of systematics and the correlations between 
redshift bins we refer to [27|]. We perform a likelihood 
analysis for the amplitude of the tensor modes At and 
marginalize over the following parameters: the cosmo- 
logical constant in the range 0.6 < fl\ < 0.9, the scalar 



spectral index 0.8 < n s < 1.3 and a Hubble parameter 
0.6 < h < 0.9. We also vary cr 8 between 0.7 and 1.1. We 
assume a spatially-flat cosmological model, and matter 
density changes as fl\ is varied with the condition that 
Q m = 1 — Q\. To obtains constrains on the tensor am- 
plitude, we only make use of large scale, 10 < I < 40, 
WMAP data. In order to match the region of angu- 
lar scales probed by the ISW-correlation data and avoid 
biasing due to lack of power in WMAP quadrupole to 
the octopole, we have excluded the region of multipolc 
I < 10. We also vary the baryon density, as a function 
of h, with a fixed value for Sl^h 2 of 0.02, but found little 
variation in our results to variations in f^; This is due 
to the fact that the analysis excludes acoustic peaks that 
are more sensitive to the exact baryon density, relative to 
matter, while at large scales, corresponding to the mul- 
tipole range above, there is little or no sensitivity in the 
anisotropy spectrum to parameter Of,. 

We plot the results in Fig which shows the 1 and 2a 
contours in the — r plane. As we may notice current 
data constrain the tensor contribution to r < 0.59 at 2a 
confidence level. Since we have only 5 weak detections of 
the ISW-galaxy cross-correlation, which may be affected 
by unknown systematics, one should be careful in a di- 
rect comparison with other constraints in the literature. 
Nonetheless the results of this analysis already illustrates 
the potential of this method which captures all the infor- 
mation one can extract on the tensor amplitude given the 
knowledge of the galaxy distribution. In fact, the previ- 
ous standard approach combines only the CMB and mat- 
ter power spectrum. Similarly we use the galaxy power 
spectrum to constrain the bias parameter, but with the 
addition of the CMB-LSS cross-correlation. The only 
way to improve this method would be to consider the 
cross-correlation as a function of galaxy luminosity and 
redshift, since this would allow a better separation of the 
uncertainties related to bias and dark energy respectively. 

As shown in Fig. 1, the constraint on the tensor-to- 
scalar ratio depends on Oa- In particular larger val- 
ues of Q\ allow a slightly greater tensor contribution. 
This trend is expected, in fact for low A values the ISW- 
correlation is smaller, therefore to match the observed 
cross-correlation larger values of the scalar amplitude are 
favoured and thus strongly limiting the tensor contribu- 
tion when combined with the CMB normalization. On 
the contrary, increasing the dark energy density causes a 
larger ISW-correlation signal which will favour low val- 
ues of A s and consequently a larger contribution of the 
tensor modes. If future dark energy measurements, such 
as from Type la supernovae, are to point towards val- 
ues of the cosmological constant around Oa ~ 0.5 — 0.6, 
the current cross-correlation data already constrain the 
GW contribution to be less than 10% of the scalar am- 
plitude. Beyond dark energy, we also found a minimal 
correlation of our constraint with the spectral index ns 
over the range of values considered. Given the limited 
information from the cross-correlation we did not vary 
other parameters or included a running to the tilt. 
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It is clear that accounting for the CMB-LSS cross- 
correlation provide better limits on the tensor-to-scalar 
ratio than those derived from WMAP data alone, from 
which r < 0.9. However the improvement does not in- 
crease significantly as one includes further information 
on the CMB-LSS cross-correlation. To understand this 
point, consider the CMB temperature fluctuation at large 
angular scales, 

AT = AT SW + AT ISW + AT GW . (6) 

Assuming that all other cosmological parameters are 
known except for the amplitude of the tensor perturba- 
tion At, the limit on the tensor-to-scalar ratio is 




where 

ac ' = V^TTT (°? W + C > SW + 2^ SW ~ ISW ) , (8) 

and we have assumed cosmic variance limited measure- 
ments. Taking the full anisotropy spectrum (which in- 
cludes contributions beyond SW and ISW effects), the 
above estimate leads to a cosmic variance limit on the 
tensor-to-scalar ratio of 0.06 at the one-sigma confidence 
level. In order to improve this limit further, the CMB- 
LSS cross-correlation can be used to remove the ISW 
contribution. The extent to which the signal improves 
depends on how well Cj sw , and to a lesser extent the 
correlation between SW and ISW effects Cf W - ISW can 
be removed from the data. If a fraction c of the ISW 
effect, in tempertature maps, can be accounted for by 
the CMB-LSS cross-correlation, then the improved noise 
estimate for the calculation is 

°' Cl = [Cf W + (1 - c) 2 Q ISW + 2(1 - c)Cf w " ISW ] 

where we have ignored that the removal in the tempera- 
ture maps also reduces the correlated part with ISW by 
a similar fraction though in reality, this may be different 
given the projection effects and the redshift dependence 
of the LSS maps used to clean the CMB maps from the 
ISW contribution. Assuming that one can safely remove 
a 50% contribution of the ISW effect, the improvement 
is at the level of 5%. Therefore removing the ISW allows 
us to push the tensor-to-scalar limit only to the level of 
0.056. 



While future constraints are not significant, the ISW 
removal allows us to include information related to 
other parameters, as the ultimate limit from tempera- 
ture anisotropics requires perfect knowledge on all cos- 
mological parameters other than the tensor amplitude. 
Thus, there are significant gains to be made with us- 
ing the CMB-LSS correlations to improve the tensor-to- 
scalar limit. This is similar to the improvements one see 
when combining CMB data with LSS data. Note that, in 
future, best limits on the tensor-to-scalar ratio will come 
from CMB polarization observations, rather than tem- 
perature anisotropies. It is well-known that the direct 
signal in CMB B-models presents the best opportunity 
to establish the tensor contribution, even in the presence 
of a confusion with a lensing effect j3J, |35|, Hy, H3, Hg ■ 



III. CONCLUSIONS 

Large scale CMB anisotropies involve a combination 
of the scalar and tensor fluctuations, we have argued 
that the scalar amplitude can be independently deter- 
mined through the CMB-galaxy cross-correlation, under 
the reasonable assumption that the large-scale cluster- 
ing of galaxies is not affected by primordial gravitational 
waves. Previous limits on the tensor amplitude have con- 
sidered only the matter power spectrum to determine the 
scalar amplitude. We have shown that an improved limit 
can be obtained by including the cross-correlation be- 
tween CMB and large scale structure surveys, in addi- 
tion to the galaxy power spectrum alone. Using recently 
measured cross-correlation data, we obtain a constraint 
r < 0.5 at 68 % confidence level on the tensor-to-scalar 
ratio. The data also allow us to exclude gravity waves at 
a level of a few percent, relative to the density field, in a 
low - Lambda dominated universe ~ 0.5). This tech- 
nique can account for the knowledge of additional cos- 
mological parameters which are required to determined 
the amplitude of the tensor modes and helps converging 
on the ultimate tensor-to-scalar ratio limit of 0.05 from 
temperature anisotropies alone (assuming all other cos- 
mological parameters are known). 
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